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Part  I:  The  structure  and  dynamics  of  self-assembled  di-lithium  pthalocyanine 
(Li2Pc) 

1. 1.  Introduction 

Di-lithium  phthalocyanine  (Li2Pc),  as  shown  in  Figure  1,  has  been  proposed  as  a  solid 
electrolyte  for  lithium-ion  batteries,  with  the  expected  advantage  of  providing  single-ion 
transport  characteristics  for  lithium  ions.1 
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Figure  1 :  Li2Pc  molecule. 
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Molecular  self-assembly  may  lead  to  the  formation  of  lithium  ion  conducting  channels, 
where  the  anion  matrix  of  di-lithium  phthalocyanine,  an  unsaturated  macrocyclic 
compound,  forms  the  channel.  When  such  self-assembled  structure  is  used  as  solid 
electrolyte  of  lithium-ion  batteries,  lithium  ion  transport  may  depend  on  the  electric  field 
gradient  established  between  the  electrodes,  instead  of  being  a  function  of  polymer 
segmental  motion  as  that  observed  for  oxygen-based  solid  polymer  electrolytes,  and 
therefore,  it  is  expected  that  the  temperature  dependence  for  lithium  ion  conduction  will 
be  minimized.  Also,  voltage  drop  across  the  electrolyte  will  be  minimized  even  at  high 
current  loads  because  of  the  single-ion  transport  characteristics.  The  performance  of 
Li2Pc  as  electrolyte  in  an  all-solid  state  electrochemical  cell  has  been  recently  reported.2 

In  view  of  the  potential  advantages  of  the  use  of  L^Pc  as  a  solid  electrolyte, 
understanding  the  lithium-ionic  transport  mechanisms  in  solid  L^Pc  becomes  essential. 
However,  little  is  known  about  its  crystalline  structure.  Other  metal  phthalocyanines  are 
known  to  be  polymorphic.3'13  Particularly,  a  great  deal  of  research  has  been  done  on  the 
polymorphism  of  lithium  phthalocyanine  because  of  its  potential  use  as  an  oximetry 
probe  and  its  interesting  electrical  and  magnetic  behavior.  ’  Depending  on  the 
preparation  conditions,  LiPc  can  crystallize  in  three  different  structures,  namely,  a-,  P-, 
and  x-forms.  The  a-  and  P-forms  of  LiPc  show  monoclinic  unit  cells  with  parameters  a  = 
2.57  nm,  b  =  0.38  nm,  c  =  2.36  nm,  and  p  =  91.0°  and  space  group  C2/c,14  and  a  =  1.94 
nm,  b  =  0.49  nm,  c=  1.40  nm,  and  p  =  120.36°  and  space  group  P2i/C,15  respectively.  The 
X-form  shows  a  tetragonal  unit  cell  with  parameters  a  =  b  =  1.385  nm  and  c  =  0.65  nm 
and  the  space  group  P4/mcc.9  These  three  polymorphs  showed  different  properties 
correlated  to  the  nature  of  molecular  packing.10 

X-Ray  spectra  of  L^Pc  were  obtained  at  Argonne  National  Laboratory16  and  at  the  Air 
Force  Laboratory  where  L^Pc  solid  state  density  was  also  determined.  We  have  used 
this  information  to  formulate  models  that,  in  conjunction  with  ab  initio  data  for  the  single 
molecule  and  dimers  are  input  for  classical  molecular  dynamics  (MD)  simulations  used  to 
investigate  self-assembled  structures  and  determine  their  corresponding  X-Ray  spectra. 
We  focused  on  the  characterization  of  structural  and  dynamical  properties  of  a  self- 
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assembly  of  LfPc  and  its  potential  to  provide  lithium-ion  conducting  channels  for 
lithium-ion  batteries.  Calculated  X-Ray  spectra  and  lithium-ion  diffusion  coefficients  are 
compared  to  experimental  results.  The  possibility  of  existence  of  ion-channels  in  the  solid 
state  structure  is  investigated  through  analyses  of  configurations  derived  from  the  MD 
simulations. 

I.  2.  Methodology:  Molecular  dynamics  simulations. 

Four  primitive  unit  cells  (designated  as  a-,  P-,  £-,  and  y-forms)  were  built  based  on 
several  sources  of  data  for  this  and  related  systems:  1)  the  analysis  of  experimental  X- 
Ray  spectra  of  crystalline  Li2Pc,16(for  the  s-form);  2)  the  experimental  density  of  solid 
Li2Pc  at  room  temperature  (for  the  a-,  P-,  and  y-forms);  3)  the  reported  polymorphism 
of  various  metal  phthalocyanines7,8,10  (for  the  a-,  P-,  and  y-forms).  The  cell  parameters  of 
the  four  configurations  along  with  all  simulation  details  have  been  reported.18  Among 
these  polymorphs,  the  a-,  P-,  and  s-form  structures  have  monoclinic  unit  cells,  and  the  y- 
form  has  a  tetragonal  unit  cell.  The  primitive  unit  cells  for  the  a-,  P-,  and  y-forms  have 
the  same  dimensions  as  the  reported  structures  for  the  polymorphs  of  LiPc  except  that  the 
shortest  cell  dimension  was  elongated  to  fit  the  experimental  density  of  LfPc  (1.40 
g/cm  ).  The  cell  parameters  of  the  s-form  were  chosen  to  match  the  crystal  d-spacing  that 
best  reproduces  the  experimental  X-Ray  data,  as  calculated  with  the  powder  indexing 
program  CRYSFIRE  2002. 19  The  DL  POLY  program,20  version  2.13,  was  used  in  ah 
simulations.  MD  simulations  were  run  in  the  microcanonical  ensemble  (NVE)  for  1200 
ps  with  a  time  step  of  0.001  ps  at  300  K.  Equilibration  runs  of  500  ps  were  performed 
before  collecting  averages  in  production  runs  of  700  ps.  During  the  equilibration  phase, 
the  interval  for  scaling  velocities  according  to  the  selected  temperature  was  set  to  0.002 
ps,  in  the  production  phase  the  scaling  is  disconnected,  but  the  temperature  remains  about 
constant.  The  cutoff  radius,  beyond  which  intermolecular  interactions  of  the  real  space 
part  of  the  long-range  electrostatic  and  the  van  der  Waals  potentials  were  set  to  zero,  was 
chosen  as  10.0  A,  which  corresponds  to  half  the  minimum  cell  length. 

I.  3.  Results  and  Discussion 

1.3.1  Structure.  The  final  structures  (after  1200  ps  total  simulation  length)  of  the  four 
forms  of  L^Pc  (Table  1)  are  shown  in  Figure  2.  The  corresponding  Li-Li  intramolecular 
distance  (for  a-,  (1-,  s-,  and  y-forms)  are  also  computed  as  averages  and  reported  in  Table 
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7.  The  MD  results  (Figure  2)  indicate  that  molecules  in  s-form  and  %-form  are  shifted 
from  each  other  in  two  adjacent  layers,  while  those  in  the  a-form  and  P-form  are 
staggered.  In  the  shifted  structure,  lithium  atoms  locate  between  the  ring  center  of  one 
molecule  and  one  of  the  benzene  rings  of  the  second  molecule,  whereas  in  the  staggered 
structure,  both  lithium  atoms  locate  between  the  ring  centers  of  the  two  molecules.  The 
separations  between  layers  have  been  computed  (from  the  last  MD  configuration  of  each 
case)  as  averages  over  measured  distances  between  all  the  equivalent  two  adjacent  layers 
within  the  simulation  cell.18  Comparison  of  these  distances  with  those  of  (Li2Pc)2 
calculated  by  B3LYP/6-31G(d)  (3.72  A  for  the  staggered  form  and  3.77  A  for  the  shifted 
form)  indicates  good  agreement  between  MD  and  DFT  regarding  the  staggered  forms, 
whereas  shorter  distances  are  found  by  MD  in  the  shifted  configurations.18 


(c)  (d) 

Figure  2:  Structures  resulting  from  MD  after  1200ps  of  production  time 


(a)  a-form,  (b)  (3-form,  (c)  8- form,  (d)  x-form 
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Another  difference  between  the  DFT  dimer  structures  and  the  MD  results  arises 
from  the  intramolecular  distance  between  two  lithium  atoms,  which  are  elongated  after 
MD  (2.2-2. 6  A  for  s-,  and  %- forms,  and  2.4-2. 9  A  for  a-,  and  P-forms),  whereas  1.99  A, 
2.23  A,  and  2.56  A  separations  were  found  after  B3LYP/6-31G(d)  geometry  optimization 
of  a  single  molecule  as  well  as  in  the  shifted  and  staggered  dimers,  respectively.  The 
longer  Li-Li  intramolecular  distances  obtained  from  the  MD  simulations  in  comparison 
with  the  DFT  results  point  to  the  influence  of  collective  effects  of  the  much  larger 
ensemble  of  molecules  in  the  MD  simulations,  but  also  they  may  be  a  reflection  of  a 
weakening  of  the  Li-N  bonds  which  would  facilitate  the  mobility  of  the  lithium  atoms  in 
the  self-assembled  structure.  On  the  other  hand,  the  small  energy  difference  found 
between  the  average  total  energies18  is  an  indication  of  possible  polymorphism  of  L^Pc. 

1.3.2  X-ray  spectra.  PowderCell  Version  1.0,  a  program  for  exploring  and  manipulating 
crystal  structures  and  calculating  X-Ray  powder  patterns  was  used  to  calculate  X-Ray 
spectra  corresponding  to  the  simulated  structures.  This  program  uses  as  input  primitive 
cell  parameters  and  structure  (atomic  positions)  and  simulates  an  X-Ray  pattern  using 
Bragg’s  law,  weighting  the  contributions  of  each  atom  with  calculated  structure  factors. 
We  used  as  input  the  coordinates  from  the  final  molecular  structure  of  the  MD 
simulations  in  the  primitive  cells  of  each  of  the  four  crystalline  forms,  along  with  their 
respective  primitive  cell  parameters,18  to  obtain  simulated  X-Ray  powder  diffractograms. 

Figures  3a  and  b  show  experimental  X-ray  diffraction  spectra  of  L^Pc 
synthesized  in  Argonne  National  Laboratory  (dried  at  220  °C  in  vacuum),  and  received 
from  Aldrich  (dried  at  160  °C  in  vacuum)  respectively.  Comparing  the  two  experimental 
spectra,  we  observe  that  in  3  a  several  sharp  peaks  appear  at  an  angle  20  around  20°,  but 
only  a  broad  peak  is  detected  in  3b,  where  some  small  peaks  disappear.  We  speculate  that 
the  experimental  X-ray  spectra  may  reflect  a  mixture  of  L^Pc  polymorphs,  some  of  them 
less  crystalline  than  others;  the  presence  of  impurities  may  also  add  new  features  to  the 
spectrum,  as  observed  in  other  metal  phthalocyanines. 13,22  In  addition,  different  methods 
to  synthesize  L^Pc  may  give  different  mixtures  of  LFPc  polymorphs.  It  is  also  possible 
that  structural  changes  may  exist  because  of  the  different  drying  temperatures. 
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Figure  3:  Comparison  of  simulated  X-ray  with  experimental  data 

(a)  Li2Pc  synthesized  in  Argonne  National  Laboratory  (dried  at  220  °C  vacuum) 

(b)  Li2Pc  from  Aldrich  (dried  at  160  °C  vacuum) 

(c)  a-form,  (d)  (3-form,  (e)  e-form,  (f)  %-form 


Figures  3c  to  f  show  the  simulated  X-ray  spectra  for  the  four  investigated 
crystalline  forms.  In  our  comparative  analysis,  we  notice  that  internal  molecular 
reorientations  of  the  Li2Pc  molecules  in  a  primitive  cell  may  cause  enhancement, 
reduction,  or  even  disappearance  of  some  of  the  peaks  that  would  arise  because  of  the 
primitive  cell  dimensions,  and  the  possible  influence  of  polymorphism  in  the 
experimental  data  may  obscure  a  direct  comparison  of  the  calculated  and  experimental 
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structures.  However,  we  perform  a  direct  comparison  of  each  of  the  calculated  spectrum 
to  the  experimental  data  (Figs.  3a  and  b),  assuming  that  the  latter  are  single  crystalline 
phases. 

A  reasonable  agreement  in  the  peak  positions  is  found  for  the  s-form  (Figure  3e), 
not  surprisingly  since  the  parameters  of  the  primitive  cell  of  the  s-form  were  selected  to 
reproduce  the  d-spacing  of  the  experimental  diffractogram  of  Figure  3  a;  however  the 
agreement  is  less  satisfactory  for  angles  20  >  20°.  The  correspondence  is  much  better  for 
the  (3-form  which  has  similar  dimensions  to  the  s-cell.  Moreover,  comparison  between 
the  relative  intensities  of  the  simulated  diffraction  peaks  with  the  experimental  ones, 
points  to  a  better  agreement  of  the  P-form  over  the  s-form.  The  experimental  spectra  in 
Figures  3a  and  b,  and  that  of  the  calculated  p-phase  (Figure  3d)  show  also  a  close 
similarity  to  the  reported  spectra  for  the  p-phase  of  H2Pc  and  MgPc. 

On  the  other  hand,  the  simulated  X-ray  spectrum  of  the  a-form  of  Li2Pc  exhibits 
few  and  rather  broad  diffraction  peaks,  compared  to  other  three  Li2Pc  simulated 
polymorphs.  This  feature  is  indicative  of  a  low  crystallinity  of  the  a-form,  which  was  also 
observed  for  the  a- forms  of  other  metal  phthalocyanines.  ’  The  fusion  of  the  two  peaks 
appearing  in  the  range  5  <  20  <  10°  is  another  characteristic  of  the  a-phase,  also  detected 
in  other  metal  and  non-metal  phthalocyanines.  ’  In  comparison  with  Figures  3a  and  b, 
we  discard  the  contribution  of  the  a-phase  to  the  experimental  spectrum. 

The  simulated  X-ray  spectrum  of  the  x-form  shows  a  well  defined  crystalline 
structure.  However,  although  some  peaks  appear  at  the  same  position  of  the  experimental 
ones,  most  of  them  do  not,  and  probably  the  x-form  is  not  the  predominant  form  that 
contributed  to  the  experimental  data.  However,  the  peak  of  highest  intensity  for  this 
structure  appears  a  value  20  of  18.8  0  in  coincidence  with  the  broad  peak  of  Figure  3b  that 
makes  the  main  difference  between  the  experimental  structures  Figures  3a  and  b.  This 
aspect  points  out  to  the  possibility  of  a  mixture  of  the  most-likely  P-form  with  the  x-form 
in  the  experimental  (Figure  3b)  spectrum.  In  summary,  the  results  in  Figure  3  indicate 
that  the  best  agreement  between  experimental  and  calculated  results  is  found  for  the  p- 
phase. 

1.3.3  Ion  channel  formation  and  ionic  diffusion.  Looking  at  the  self-assembly  molecular 
organization,  an  alignment  of  lithium  ions  is  found  that  can  be  associated  with  the 


7 


expected  lithium  ion  conducting  channels.  We  can  visualize  these  conducting  channels  in 
a-form,  P-form,  and  s-form  when  combining  the  top  (not  shown)  and  side  views  of  the 
unit  cells,18  or  rotating  the  unit  cell  of  x-form.  These  chains  of  ions  are  oriented  along  the 
same  direction  across  layers. 


To  investigate  lithium  mobility,  we  calculated  the  velocity  autocorrelation 
function  (VAF),  which  is  related  to  the  macroscopic,  phenomenological,  self-diffusion 
coefficient  D  through  the  Green-Kubo  formula,24  written  as  the  time  integral  of  a 
microscopic  time-correlation  function: 


00 

D  =  \{v(t).v(t0))  dt 
0 


(1) 


where  the  function  in  brackets  is  a  measure  of  the  projection  of  the  particle  velocity  at 
time  t  onto  its  initial  value,  averaged  over  all  initial  conditions  to.  The  velocity 
autocorrelation  function  of  lithium  ions  shows  oscillations  in  the  time  domain,  indicating 
the  rattling  motion  of  lithium  ions  in  the  “cage”  of  their  nearest  neighbors.18  The 
corresponding  conductivity  of  lithium  ions  was  calculated  via  the  Nemst-Einstein 
equation, 
cDz2F2 


K  = 


RT 


(2) 


where  the  c  is  the  concentration  of  lithium  ions,  D  is  the  ionic  diffusion  coefficient,  z  is 
the  charge  on  lithium  ions,  F  is  the  Faraday  constant,  R  is  the  gas  constant,  and  T  is  the 
absolute  temperature. 


Table  1:  Calculated  diffusion  coefficient  and  conductivity  of  lithium  ions  in  various  Li2Pc  crystalline 

forms  at  300K. 


Li2Pc  crystalline  form 

Diffusion  coefficient 
(cm2/s) 

Conductivity* 

(S/cm) 

a 

5.94xl0“9 

2.65x10'5 

p 

1.35xl0-8 

6.03xl0"5 

8 

2.04xl0'8 

9.79xl0'5 

X 

1.88xl0'8 

8. 40x1  O'5 

*Experimental  ionic  conductivities  of  5.1  to  8.9  xlO"4  S/cm  were  obtained  at  300  K.17 
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The  results  shown  in  Table  1  illustrate  the  differences  in  ionic  mobility  imposed  by  the 
various  crystal  structures,  with  the  less  crystalline  a-form  having  the  smallest  self- 
diffusion  coefficient  and  the  s-form  the  highest.  A  transport  mechanism  for  lithium  ions 
in  Li2Pc  has  been  suggested,  whereby  instead  of  depending  on  polymer  segmental  motion 
as  in  the  PEO  system,  the  lithium  ion  diffusion  would  depend  mainly  on  the  electric 
field  gradient  in  the  system.  Values  of  diffusion  of  lithium  ions  in  PEO/LiCICE  mixtures 
were  reported  to  be  5.2  x  10"9  cm2/s  at  349  K,26  which  are  below  our  calculated  values  at 
300  K.  Compared  to  the  experimental  ionic  conductivity  in  L^Pc  (Table  1),  the 
simulated  lithium  ion  conductivity  results  are  about  one  order  of  magnitude  smaller.  This 
difference  is  explained  considering  that  in  the  simulation  ions  move  via  self-diffusion 
determined  by  the  field  created  by  the  self-assembly  of  L^Pc  molecules,  while  in  the 
experiment  diffusion  across  the  electro lyte/electrode  interface  is  driven  by  a  difference  of 
chemical  potential  and  by  an  electric  field  gradient  established  between  the  electrodes  of 
an  electrochemical  cell.  Work  to  incorporate  the  effect  of  an  electric  field  in  the 
simulations  is  in  progress,  and  will  be  reported  elsewhere.  Further,  we  note  that  equation 
(1)  has  the  implicit  assumption  that  ions  at  low  concentrations  do  not  interact  with  each 
other  and  therefore  their  motion  is  uncorrelated.  However,  this  is  an  oversimplification 
of  this  system  that  has  a  strong  long-range  Coulombic  component.  A  more  accurate 

97 

equation  would  relate  the  conductivity  to  the  charge  flux  autocorrelation  function, 
where  the  charge  flux  is  calculated  by  the  sum  of  charge  times  velocity  for  each  ion. 
Alternatively,  one  could  calculate  Maxwell- Stefan  diffusion  coefficients,  where  the  mass 
fluxes  are  linearly  related  to  the  gradient  of  chemical  potential  instead  of  being 

98 

proportional  to  the  gradient  of  concentration. 

Part  II:  Properties  of  self-assembled  LEPc  at  the  LhPc-pyrite  interface 
II.  1.  Background 

The  crystalline  structure  and  self-assembly  characteristics  of  di-lithium  phthalocyanine 
were  investigated  in  the  first  part  of  our  work.  However,  since  lithium  ions  diffuse  from 
the  electrolyte  to  the  electrode,  a  good  interfacial  contact  is  essential.  For  pyrite 
cathodes,  a  thin  film  of  L^Pc  is  grown  over  the  surface  of  pyrite.  Pyrite  has  a  simple 
cubic  structure  with  space  group  Pa3  and  a  unit  cell  with  length  parameter  of  5.4281  A 
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containing  four  FeS2  units  define  this  system.  Figures  4a  and  b  show  the  bulk  pyrite  and 
pyrite  (100)  surface.  In  bulk  pyrite,  each  Fe  is  coordinated  by  6  S  atoms  and  each  S  by  3 
Fe  atoms  and  1  S  atom.  In  the  (100)  surface,  Fe  is  coordinated  by  5  S  atoms;  S  by  2  Fe 
atoms  and  1  S  atom. 


(a) 


Figure  4:  Structural  models  of  bulk  pyrite  (a),  and  pyrite  (100)  surface  (b).  The  yellow  spheres  are  the 

S  atoms,  and  the  red  ones  are  the  Fe  atoms. 


We  use  molecular  dynamics  simulations  to  analyze  the  pyrite/Li2Pc  solid-solid  interface, 
with  the  aim  of  obtaining  structural  details  of  a  self-assembled  Li2Pc  electrolyte  film  in 
contact  with  a  pyrite  cathode  surface.  Also,  we  try  to  understand  the  solid-solid 
interfacial  contact  between  Li2Pc  and  the  pyrite  surface  and  investigate  the  effect  of  the 
presence  of  the  pyrite  cathode  on  the  formation  of  the  lithium-ion  conduction  channels  in 
the  Li2Pc  phase.  By  investigating  the  interfacial  structure  of  the  system,  several 
structural  and  dynamical  properties  are  analyzed,  including  the  radial  distribution 
functions,  and  adsorption  of  Li2Pc  on  the  pyrite  (100)  surface. 

II.  2.  Simulation  details 

The  simulation  system  consists  of  a  5-layer  slab  of  a  pyrite  exposing  a  (100)  surface 
in  contact  with  a  mono-layer,  a  bi-layer,  or  multi-layers  of  Li2Pc  located  in  a  tetragonal 
unit  cell,  which  is  subjected  to  periodic  boundary  conditions.  The  5-layer  slab  of  pyrite 
contains  250  FeS2.  During  the  simulation,  the  pyrite  atoms  in  the  slab  are  fixed.  The 
atomic  charges  of  +0.64  for  Fe,  and  -0.32  for  S  atoms  were  taken  from  published  results. 
The  van  der  Waals  force  field  parameters  were  taken  from  the  literature.29  The  simulation 
cell  has  dimensions  of  50.00  A  x  27.14  A  x  27.14  A. 
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So  far,  the  initial  configuration  of  Li2Pc  in  contact  with  pyrite  (100)  surface  is  of  the  x- 
form  previously  proposed.18  Two  initial  orientations  of  the  x-form  Li2Pc  were  studied, 
where  the  stacking  axis  of  the  self-assembled  Li2Pc  is  either  perpendicular  or  parallel  to 
the  pyrite  (100)  surface  respectively.  For  initial  configurations  with  stacking  axis 
perpendicular  to  the  surface,  5  cases  with  systems  containing  1,  2,  4,  8,  and  16  Li2Pc 

were  simulated.  The  LFPc  molecules  are  arranged  in  layers 
parallel  to  the  surface,  with  a  maximum  number  of  four  Li2Pc 
molecules  in  each  layer.  Figure  5  shows  the  initial  configuration 
of  the  pyrite  slab  with  4  Li2Pc  molecules  in  contact  with  the 
surface,  with  an  initial  separation  of  2.5  A  between  the  Pc  ring 
and  the  surface. 

For  the  cases  containing  one  and  two  Li2Pc  molecules  the 
separation  distance  from  the  surface  is  kept  the  same.  For  the 
cases  contain  8  and  16  Li2Pc  molecules,  additional  layer  (s)  were 
added  to  the  top  of  the  first  layer  shown  in  Figure  5,  with 
separation  of  3.25  A  between  Pc  rings.  For  the  simulations 
where  the  stacking  axis  of  the  self-assembled  Li2Pc  was  parallel 
to  the  pyrite  (100)  surface,  four  cases  with  systems  containing  1, 
2,  8,  and  16  Li2Pc  molecules  were  simulated.  Figure  6a  shows 
the  top  view  of  the  initial  configuration  of  the  pyrite  slab  in  contact  with  a  self-assembly 
of  8  Li2Pc  molecules  covering  one  half  of  the  pyrite  (100)  surface.  Figures  6  (b),  (c),  and 
(d)  show  the  systems  containing  1,2,  and  16  Li2Pc. 

The  DL  POLY  program,  version  2.13,  was  used  in  all  MD  simulations.  MD 
simulations  were  run  in  the  microcanonical  ensemble  (NVE)  for  800  ps  with  a  time  step 
of  0.001  ps.  Equilibration  runs  of  300  ps  were  performed  before  collecting  averages  in 
production  runs  of  500  ps.  During  the  equilibration  phase,  the  interval  for  scaling  of 
velocities  according  to  the  selected  temperature  was  set  to  0.002  ps;  in  the  production 
phase  the  scaling  is  disconnected,  but  the  temperature  remains  about  constant.  The  cutoff 
radius,  beyond  which  intermolecular  interactions  of  the  real  space  part  of  the  long-range 
electrostatic  and  the  Van  der  Waals  potentials  were  set  to  zero,  was  chosen  as  10.0  A. 


Figure  5.  Initial 
configuration  of  the 
pyrite  slab  with  4 
Li2Pc  molecules  (a 
monolayer)  in  contact 
with  the  pyrite  (100) 
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Figure  6:  the  initial  configuration  ofpyrite  slab  with  stacking  axis  of  Li2Pc  parallel  to  (100)  surface: 
(a)  8  Li2Pc,  (b)  1  Li2Pc,  (c)  2  Li2Pc,  (d)  16  Li2Pc 


II.  3.  Results  and  discussion 

II.  3.1  Initial  stacking  axis  of  Li2Pc  self-assembly  is  perpendicular  to  the  pyrite  (100) 
surface 

Figure  7  shows  the  side  view  of  the  final  configurations  of  systems  with  various  film 
thicknesses  of  LCPc.  Compared  to  the  initial  separation  distance  between  the  Pc  ring  and 
the  (100)  surface,  the  mean  separation  distance  increased  from  2.5  A  to  3.5  A  after  the 
system  equilibrates.  This  is  because  of  the  main  interactions  between  pyrite  and  LCPc, 
i.e.,  the  electrostatic  attraction  between  Fe  and  N,  Li  and  S,  Li  and  N,  as  well  as  the 
repulsion  between  Fe-Li,  S-N,  and  Li-Li.  For  the  systems  containing  more  than  one  L^Pc 
molecule,  shifted  dimers  were  formed. 
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Figure  7:  Final  configuration  of  systems  with  stacking  axis  of  Li2Pc  perpendicular  to  the  pyrite  (100) 
surface  with  varying  number  of  molecules  in  the  simulation  cells:  (a)  1  Li2Pc,  (b)  2  Li2Pc,  (c)  4  Li2Pc, 

(d)  8  Li2Pc,  (e)  16  Li2Pc 


Figure  7  (b)  shows  two  neighbor  cells,  illustrating  the  formation  of  the  shifted  dimer. 
Clearly,  we  can  observe  that  the  system  with  higher  Li2Pc  density  retains  the  %-form  of 
the  Li2Pc  crystalline  structure.  From  the  top  view  in  Figures  8a  to  d,  we  can  see  that  the 
lithium  ions  which  are  closest  to  the  pyrite  (100)  surface  are  adsorbed  on  sites  close  to  the 
top  of  S  atoms.  It  is  understandable  that  lithium  ions  are  not  exactly  on  the  top  of  S 
atoms,  considering  the  site  of  Li2Pc  molecule  and  thus  the  overall  repulsion  and  attraction 
interactions  between  different  ion  pairs  mentioned  previously.  For  example,  the  location 
of  Li-ions  are  mainly  decided  by  the  interactions  of  Li-Li,  Li-Fe,  Li-S,  and  Li-N.  As  the 
number  of  Li2Pc  molecules  increases  in  the  system  up  to  16  (where  the  film  is  composed 
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of  four  layers),  we  begin  to  see  the  existence  of  potential  lithium  ion  conducting  channels 
that  were  observed  in  our  previous  work  of  bulk  self-assembled  Li2Pc.  Figure  8  (e) 
shows  such  potential  conducting  channels  in  four  unit  cells. 


(d)  (e) 

Figure  8  (a)  to  (d):  Top  view  of  final  configurations  of  the  system  with  1  to  8  Li2Pc  molecules;  (e) 
view  of  potential  lithium  ion  conducting  channels  of  the  system  with  16  Li2Pc  molecules.  For  better 

visualization,  lithium  ions  are  highlighted. 

Figures  9  (a)  and  (b)  show  the  radial  distribution  functions  of  Li-S  and  Fe-N.  A  sharp 
first  peak  in  Figure  9  (a)  indicates  a  strong  interaction  between  the  Li  ions  and  the  S 
atoms  from  the  surface.  The  average  shortest  distance  is  2.65  A.  The  decrease  of  the  first 
peak  magnitude  is  due  to  the  relatively  smaller  atomic  fraction  of  Li  atoms  on  the  surface 
as  the  Li2Pc  density  increases.  The  position  of  the  first  peak  shifts  to  the  left  as  the 
density  of  Li2Pc  increases,  indicating  enhanced  S-Li  interactions.  The  negatively  charged 
N  atoms  of  Li2Pc  interact  with  the  positively  charged  Fe  atoms  of  the  substrate;  however, 
the  shortest  distance  Fe-N  is  larger  than  4  A,  showing  that  the  molecule  attaches  to  the 
surface  via  the  Li-S  interaction. 
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Figure  9:  Radial  distribution  function  of  initial  Li2Pc  stacking  axis  perpendicular  to  pyrite  (100) 

surface,  (a)  Li-S,  (b)  Fe-N 


In  order  to  test  the  formation  of  shifted  Li 2  Pc  dimers  above  the  pyrite  surface  is  not 
caused  by  the  small  pyrite  (100)  surface  contacted  with  Li2Pc,  we  used  another  5  layer 
pyrite  slab  with  (FeS2)64o  and  4  Li2Pc.  The  pyrite  (100)  surface  has  dimensions  of  43.42 
A  x  43.42  A.  Initially,  four  Li2Pc  were  located  on  the  same  plane  parallel  to  the  pyrite 

(100)  surface,  with  a  separation  distance  of  2.5  A 
between  each  Pc  ring  and  the  surface.  Compared  to 
previous  systems,  the  surface  has  a  much  larger  area  than 
that  occupied  by  the  four  Li2Pc  molecules.  Figure  10 
shows  the  final  configuration  of  the  system.  The  final 
configuration  indicated  that  one  of  the  Li2Pc  molecules 
moved  away  from  the  surface  and  formed  a  second  layer. 


By  checking  the  distance  from  the  Li  ions  belonging  to  two  Li2Pc  molecules  in  two 
different  layers,  it  was  found  that  one  pair  had  the  value  of  5.28  A,  which  corresponded 
to  the  intermolecular  Li-Li  distance  of  the  shifted  dimer.  Other  two  Li-Li  pairs  were  more 
than  10  A  away.  The  result  shows  that  the  formation  of  shifted  dimer  is  favorable  and 

stable.  The  formation  of  the  shifted  dimer  is  energetically  favorable  according  to  DFT 

1  8 

calculations  as  well. 


Figure  10:  (FeS2)64o  with  4  Li2Pc 
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11.3.2  Initial  stacking  axis  ofLi2Pc  parallel  to  the  pyrite  (100)  surface 

Recent  studies  have  found  that  most  dilithium  octacyanophthalocyanine 
[Li2Pc(CN)s]  molecules  in  the  Langmuir-Blodgett  film  were  stacked  in  a  face-to-face 
aggregated  state  with  stacking  axis  parallel  to  a  ITO  surface.  Figures  1 1  (a)  through  (d) 
show  the  final  configurations  of  the  systems  with  initial  stacking  axis  of  Li2Pc  parallel  to 
the  pyrite  (100)  surface,  with  1,  2,  8,  and  16  Li2Pc  molecules  respectively. 
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Figure  11:  Final  configuration  of  systems  with  initial  stacking  axis  of  Li2Pc  parallel  to  the  pyrite  (100) 
surface,  (a)  1  Li2Pc,  (b)  2  Li2Pc,  (c)  8  Li2Pc,  (d)  16  Li2Pc 


For  the  systems  containing  one  or  two  Li 2 Pc,  the  final  configurations  show  that  the  whole 
ring  of  Li2Pc  rotated  and  the  stacking  axis  of  the  self-assembled  system  changed  from  the 
parallel  to  the  perpendicular  position.  This  is  due  to  the  strong  interaction  between  Li  and 
S  atoms,  as  well  as  to  the  available  free  space  for  Li2Pc  molecules  to  rotate.  This,  on  the 
other  hand,  indicates  that  the  configuration  with  stacking  axis  of  Li2Pc  perpendicular  to 
the  pyrite  (100)  surface  is  more  stable.  Also,  we  clearly  observed  the  formation  of  a 
shifted  Li2Pc  dimer  for  the  system  with  two  Li2Pc  molecules.  However,  for  the  systems 
containing  8  and  16  Li2Pc,  we  did  not  observe  the  rotation  of  the  stacking  axis  from 
parallel  position  to  perpendicular  position.  Li2Pc  molecules  in  the  system  with  8  Li2Pc 
did  rotate  a  little  bit  and  the  Pc  rings  were  no  longer  completely  perpendicular  to  the 
pyrite  surface,  but  with  an  edge-on  configuration  on  the  pyrite  surface.  The  configuration 
of  Li2Pc  in  this  system  is  similar  to  the  (3-form  crystalline  structure  we  proposed 
previously.18  For  the  system  containing  sixteen  Li2Pc  molecules,  the  Pc  rings  of  the  final 
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configuration  were  relatively  perpendicular  to  the  pyrite  surface.  Although  the  strong 
interaction  between  Li  ions  and  S  ions  would  make  the  rotation  of  Li2Pc,  the  unavailable 
free  space  was  the  main  reason  that  Li2Pc  could  not  rotate  since  the  size  of  Li2Pc  is 
relatively  large.  The  final  configuration  of  Li2Pc  in  this  system  was  the  x-form  crystalline 
structure. 

As  we  mentioned  above  that  the  final  configurations  of  systems  with  1  and  2  Li2Pc 
molecules  were  very  similar  to  the  systems  with  the  systems  of  initial  stacking  axis  of 
Li2Pc  perpendicular  to  pyrite  surface.  Figures  12  (a)  and  (b)  show  the  radial  distribution 
functions  of  Li-S  and  Fe-N,  respectively. 


Figure  12:  Radial  distribution  function  of  initial  Li2Pc  stacking  axis  parallel  to  pyrite  (100)  surface,  (a) 

Li-S,  (b)  Fe-N 

From  Figure  12  (a),  we  can  see  that  the  average  position  of  the  Li-ions  with  respect  to  the 
S  atoms  on  the  surface  remains  equal  to  2.65  A  for  the  cases  with  1  and  2  molecules.  At 
higher  densities  of  Li2Pc,  the  number  of  molecules  in  contact  with  the  surface  is 
significantly  reduced,  revealing  a  clear  change  in  the  interfacial  configuration  of  the 
Li2Pc  assembly.  The  N  atoms  interact  with  the  Fe  atoms  with  shortest  distances  found  at 
about  4.5  A  in  Figure  12  (b),  again  we  observe  the  change  in  configuration  as  the  Li2Pc 
density  increases. 

II.  3. 3  Comparison  of  the  total  energy 

Table  2  shows  the  comparison  of  the  total  energy  for  the  various  systems.  At  constant 
Li2Pc  density,  the  energies  of  the  various  final  configurations  are  only  slightly  favorable 
to  the  case  where  the  stacking  axis  is  perpendicular  to  the  surface,  probably  because  it 
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maximizes  the  Li-S  and  Fe-N  interactions.  But  the  small  energy  difference  suggests  the 
possible  existence  of  polymorphic  structures. 


Table  2:  Comparison  of  total  energy  (kcal/mol  unit  cell) 


Number  of  Li2Pc 

Initial  stacking  axis  position  to  surface  | 

Perpendicular 

Parallel 

1 

-8.0796xl03 

-8.0648xl03 

2 

-8.4809x1 03 

-8.5448x1 03 

8 

-1.1038xl04 

-1.1059xl04 

16 

-1.4608xl04 

-1.4531xl04 

III.  Conclusions 

1  O 

Ab  initio  and  molecular  dynamics  calculations  showed  the  viability  of  using  a  large 
array  of  molecules  to  form  a  lithium  ion  conducting  channel  via  molecular  self-assembly. 
The  elongation  of  the  Li-Li  intramolecular  distances  based  on  the  molecular  dynamics 
calculations  illustrate  collective  effects  provided  by  the  self-assembled  structure  and  their 
importance  for  facilitating  fast  ionic  transport  within  the  channel. 

The  close  values  of  the  calculated  binding  energies  (B3LYP/6-31G(d))  for  the 
staggered  and  shifted  dimer  structures  suggest  the  ease  of  formation  of  various 
polymorphs  simultaneously  during  material  preparation,  although  from  the  comparison  of 
the  simulated  and  experimental  X-Ray  data  we  conclude  that  Li2Pc  is  in  a  predominant  p- 
phase.  On  the  other  hand,  calculated  ionic  conductivities  are  one  order  of  magnitude 
below  the  experimental  values,  which  is  primarily  attributed  to  the  absence  of  an  external 
driving  force  for  ionic  diffusion  in  the  simulated  self-assembled  structures. 

Molecular  dynamics  simulations  indicate  the  possibility  of  existence  of  possible 
polymorphic  structures  of  Li2Pc  over  the  pyrite  (100)  surface.  Energetically,  structures 
with  the  stacking  axis  perpendicular  to  the  surface  are  favorable,  although  the  difference 
is  not  significant.  In  the  cases  of  final  configurations  of  LTPc  with  a  stacking  axis 
perpendicular  to  the  pyrite  (100)  surface,  Li  ions  are  adsorbed  over  the  S  ions  on  the 
pyrite  (100)  surface.  The  interaction  between  Li  ions  and  S  ions  increase  as  the  number  of 
the  Li2Pc  over  the  surface  increases.  At  low  density  of  Li2Pc,  shifted  dimers  are  favorable 
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and  more  stable.  As  the  density  increases,  both  p-  and  y-form  structures  are  found  in 
different  systems  with  8  L^Pc.  At  high  LLPc  density,  only  the  x-form  structure  is  found. 
Also,  the  potential  lithium  ion  conducting  channels  can  be  realized  at  high  density. 
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